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Abstract. Density functional theory (DFT) became a universal approach to compute ground- 
state and excited configurations of many-electron systems held together by an external one-body 
potential in condensed-matter, atomic, and molecular physics. At present, the DFT strategy is 
also intensely studied and applied in the area of nuclear structure. The nuclear DFT, a natural 
extension of the self- consistent mean- field theory, is a tool of choice for computations of ground- 
state properties and low-lying excitations of medium-mass and heavy nuclei. Over the past 
thirty-odd years, a lot of experience was accumulated in implementing, adjusting, and using the 
density-functional methods in nuclei. This research direction is still extremely actively pursued. 
In particular, current developments concentrate on (i) attempts to improve the performance 
and precision delivered by the nuclear density-functional methods, (ii) derivations of density 
functionals from first principles rooted in the low-energy chromodynamics and effective theories, 
and (iii) including effects of low-energy correlations and symmetry restoration. In this study, 
we present an overview of recent results and achievements gained in nuclear density-functional 
methods. 



1. Introduction 

The Density Functional Theory (DFT) was introduced in atomic physics through the Hohenberg- 
Kohn [lj and Kohn-Sham [2] theorems. Its quantum-mechanical foundation relies on a simple 
variational concept that uses observables as variational parameters. Namely, for any Hamiltonian 
H and observable one can formulate the constraint variational problem, 

S(H - XQ) = 0, (1) 

whereby the total energy of the system E = (H) becomes a function of the observable Q = (Q), 
that is E = E(Q), provided the Lagrange multiplier A can be eliminated from functions E(X) and 
Q(X) that are obtained from the variation in Eq. (pQ) performed at fixed A. This can be understood 
in terms of a two-step variational procedure. First, Eq. (pQ) ensures that the total energy is 
minimized at fixed Q, and second, the minimization of E(Q) in function of Q gives obviously 
the exact ground-state energy Eq = ming E(Q) and the exact value Qo of the observable Q 
calculated for the ground-state wave function. 

Function E{Q) is thus the simplest model of the density functional. However, the idea 
of two-step variational procedure can be applied to an arbitrary observable or a set of 
observables, and hence the total energy can become a function E(Qk) of several observables 



Q/c, S(H — ^Zk^kQk) — =>* E — E(Qk), or a functional E[Q(q)] of a continuous set of 
observables Q(q), S(H - J dq X(q)Q(q)) = 0^E = E[Q(qj\. 

When these ideas are applied to the observable p(r) = Ya=i ~~ Y i)-> which is the local 
density of a many-body system at point r, we obtain the original local DFT, 

S(H- yW(r)p(r))=0 => E = E[p(r)], (2) 

whereby the local external potential V(r) plays the role of the Lagrange multiplier that selects a 
given density profile p(r). By the same token, the nonlocal DFT is obtained by using a nonlocal 
external potential V(r,r'), 

S(H - J dr J dr 7 y(r^)p(r,r 7 )) =0 =^ £ = #[p(r, r ')]. (3) 

In each of these cases, by minimizing the functional E[p(r)] or E[p(r, r 7 )], we obtain the exact 
ground-state energy of the many-body system along with its exact local p(r) or nonlocal p(r, r 7 ) 
one-body density. 

It is thus obvious that the idea of two-step variational principle, which is at the heart of DFT, 
does not give us any hint on which observable has to be picked as the variational parameter. 
Moreover, the exact derivation of the density functional is entirely impractical, because it involves 
solving exactly the variational problem that is equivalent to finding the exact ground state. If 
we were capable of doing that, no DFT would have been further required. Nevertheless, the 
exact arguments presented above can serve us as a justification of modelling the ground-state 
properties of many-body systems by DFT, which, however, must be rather guided by physical 
intuition, general theoretical arguments, experiment, and exact calculations for simple systems. 



2. Fundamentals 

An example that highlights connection with ab-initio theory is shown in Fig. [TJ Ref. [3], where 
binding energies calculated exactly for neutron drops confined within the Woods-Saxon potential 
of two different depths, Vq = —25 and — 35.5MeV, are compared with those corresponding to 
several Skyrme functionals [3]. Calculations of this type, performed at different depths, surface 
thicknesses, and deformations of the confining potential may allow for a better determination of 
Skyrme-functional parameters. 

In practice, the exact density functionals (Ej) or (|3j) are 
modelled as integrals of energy densities and thus 
they are called energy-density functionals (EDFs). 
They can be local functions of local densities (j4|), 
quasilocal functions of local higher-order densities (Ej), 
nonlocal functions of local densities (|6j), or nonlocal 
functions of nonlocal densities (H). 
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Figure 1. Energies of neutron drops 
confined by the Woods-Saxon potential, 
see the text. From Ref. |3J; picture 
courtesy of R.B. Wiringa. 



E\p(T)]=jdrH(p(r)), 
£[p(r)]=|drH(p(r),r(r),Ap(r), 
E[p{v)\=jdvjdv'n{p{v),p{v')), 
E[p{r,r')] =JdrJdr'n(p(r,r')). 
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Two major classes of approach that are currently used and developed in nuclear structure 
physics are based on relativistic and nonrelativistic EDFs [5JHJE]. The nonrelativistic EDFs are 
most often built as: 



Vt 



i0 



-KT(r,r') 



P*(r,r>t( r » 



lO 



(8) 



where V^ ir and V^ xc denote the EDF-generating (pseudo) potentials in the direct and exchange 
channels, respectively, and we sum over the spin-rank v = and 1 (scalar and vector) and isospin- 
rank t = and 1 (isoscalar and isovector) spherical-tensor densities [7J [8] coupled to the total 
isoscalar (superscript 0) and scalar (subscript 0) term. For example, finite-range momentum- 
independent central potentials generate the Gogny [9] or M3Y [10] nonlocal functionals (j8j) 
and zero-range momentum-dependent pseudopotentials generate the Skyrme [I] or BCP [TT] 
quasilocal functionals (EJ). 

Expression (jHj) derives from the Hartree-Fock formula for the average energy of a Slater 
determinant. However, the EDF-generating pseudopotentials should not be confused with 
the nucleon-nucleon (NN) bare or effective interaction or Brueckner G-matrix. Indeed, their 
characteristic features are different - they neither are meant to describe the NN scattering 
properties, as the bare NN force is, nor are meant to be used in a restricted phase space, as 
the effective interaction is, nor depend on energy, as the G-matrix does. Moreover, to ensure 
correct saturation properties, the EDF-generating pseudopotentials must themselves depend on 
the density. But most importantly, the generated EDFs are modelled so as to describe the exact 
binding energies and not those in the Hartree-Fock approximation, which otherwise would have 
required adding higher-order corrections based on the many-body perturbation theory. 

Only very recently, it has been demonstrated [12j [8] that the nuclear nonlocal EDFs, based 
on sufficiently short-range EDF-generating pseudopotentials, are equivalent to quasilocal EDFs. 
In Fig. [2] are compared the proton RMS radii and binding energies of doubly magic nuclei, 
determined by using the Gogny D1S EDF [13J and second-order Skyrme-like EDF SISb [12] 
derived therefrom by using the Negele-Vautherin (NV) density-matrix expansion (DME) [H]. 
One can see that already at second order, the DME gives excellent precision of the order of 1%. 
In Ref. [15], similar conclusions were also reached when comparing the nonlocal and quasilocal 
relativistic EDFs, see Fig. [3l 
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Figure 2. Proton RMS radii 
(left) and binding energies (right) 
of doubly magic nuclei, determined 
by using the Gogny D1S EDF [13] 
(squares) and second-order Skyrme- 
like EDF [12J derived therefrom 
(circles). The insets show results 
in expanded scales and scaled by 
particle numbers A. 



Figs. [Hand [5] show convergence of the direct and exchange interaction energies, respectively, 
when the Taylor and damped Taylor (DT) DMEs are performed up to sixth order [8]. The four 
panels of Fig. [5] show results obtained in the four spin-isospin channels labeled by V vt . Results of 
the DT DME [8] are compared with those corresponding to the NV [2] and PSA [16] expansions. 
It is extremely gratifying to see that in each higher order the precision increases by a large factor, 
which is characteristic to a rapid power-law convergence. The success and convergence of the 
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Figure 3. The relative deviations (in percentages) 
between the experimental and theoretical charge radii 
(upper panel) and binding energies (lower panel) of 12 
spherical nuclei, calculated with the meson-exchange 
interaction DD-ME2 and the four point-coupling 
parameter sets. In all cases, the relative deviations are 
below 1%. Starting from the meson-exchange density- 
dependent interaction DD-ME2, the equivalent point- 
coupling parametrization for the effective Lagrangian 
was derived by incorporating the density dependence of 
the parameters. The parameters of the point-coupling 
model were adjusted to reproduce the nuclear matter 
equation of state obtained with the DD-ME2 interaction 
and surface thickness and surface energy of semi-infinite 
nuclear matter. Reprinted figure with permission from 
Ref. [15J. Copyright 2008 by the American Physical 
Society. 
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Figure 4. Precision of the Taylor DME of 
the direct interaction energies calculated for 
the nonlocal Gogny D1S EDF [Hj. The nine 
nuclei used for the test are 4 He, 16 , 40,48 Ca, 
56 ' 78 Ni, 100 ' 132 Sn, and 208 Pb. From Ref. [8]. 
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Figure 5. The RMS deviations between 
the exact and approximate exchange energies 
calculated for the nine nuclei listed in the 
caption of Fig. HJ see the text. From Ref. [8]. 



DME expansions relies on the fact that the finite-range nuclear effective interactions are very 
short-range as compared to the spatial variations of nuclear densities. The quasilocal (gradient) 
expansion in nuclei works! 

3. Applications 

When compared to the experimental binding energies, the quasilocal Skyrme functional HFB-17 
[18J (Fig. E]) gives results, which have the quality very similar to those given by the nonlocal Gogny 
functional DIM [TH] (Fig. [7]). In both cases, the functional were augmented by terms responsible 
for the pairing correlations and all parameters were adjusted specifically to binding energies. 
Moreover, in both cases, by using either the 5D collective Hamiltonian approach or configuration 
mixing, theoretical binding energies were corrected for collective quadrupole correlations. The 
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Figure 6. Differences between measured |20j 
and HFB-17 [18J masses, as a function of the 
neutron number N. Reprinted figure with 
permission from Ref. |18j . Copyright 2009 by 
the American Physical Society. 



20 40 60 80 100 120 140 160 

N 

Figure 7. Differences between measured |20j 
and DIM [19] masses, as a function of the 
neutron number N. Reprinted figure with 
permission from Ref. |19| . Copyright 2009 by 
the American Physical Society. 



results are truly impressive, with the RMS deviations calculated for 2149 masses being as small 
as 798 and 581 keV, respectively. 

The problem of treating collective correlations and excitations within the DFT or EDF 
approaches is one of the most important issues currently studied in applications to nuclear 
systems. The question of whether one can describe these effects by using the functional only is 
not yet resolved. In practice, relatively simple functionals that are currently in use, require adding 
low-energy correlation effects explicitly. This can be done by reverting from the description in 
terms of one-body densities back to the wave functions of mean-field states. For example, for 
quadrupole correlations this amounts to using the following configuration-mixing states, 

\*nz,jm) = I d(3d 7 ^f K (f3^)P N P z Pj MK \*(f3^)), (9) 

J K 

where P/v, Pz, and Pjmk are projection operators on good neutron number N, proton number 
and angular momentum J with laboratory and intrinsic projections M and K. The intrinsic 

mean-field wave functions |\I/(/3,7)) correspond to one-body densities constrained to quadrupole 

deformations /3 and 7 [2T] . 

To determine variationally the mixing amplitudes /k(/3,7), one has to generalize the energy 

densities, such as those shown in Eqs. (H])-([7]), to transition energy densities that enable us to 

compute Hamiltonian kernels. For mean-field states, this can be rigorously done by using the 

Wick theorem, whereby the average energy (^\H\ty) generalizes to matrix element (\I/i|i?|\I/2) 

as [21]: 

<*|£|*) ~ jdvj dr'W(p(r,r')) for p(r,r') = (10) 
<*i|#l*2> ^ /dr|dr^( Pl2 (r,rO) f OT p 12 (r,0 = ^ (11) 

Although for densities of correlated states, which are employed in DFT or EDF methods, 
this prescription cannot be properly justified, it has been successfully used in many practical 
applications. However, even this simple prescription creates problems [22], which may require 
implementing more complicated schemes [23]. 




Figure 8. Scatter plot comparing 
the theoretical and experimental 2f 
excitation energies of the 359 nuclei 
included in the survey of Ref. |24j; 
picture courtesy of M. Bender. 



Figure 9. Theoretical 2\ excitation energies of 537 
even- even nuclei as a function of their experimental 
values [25]. Reprinted figure with permission from 
Ref. [13J. Copyright 2010 by the American Physical 
Society. 



An example of the results is shown in Fig. El where calculated 2\ excitation energies |24j 
are compared with experimental data. One obtains fairly good description for nuclei across the 
nuclear chart. Calculations slightly overestimate the data, which is most probably related to the 
fact that in this study the nonrotating mean-field states were used, see Ref. [26] and references 
cited therein. As shown in Fig. [9l this deficiency disappears when the moments of inertia of the 
5D collective Hamiltonian are determined by using infinitesimal rotational frequencies [13]. At 
present, calculations using in light nuclei the triaxial projected states of Eq. (|9j) are becoming 
possible for the relativistic (Fig. [TO]) and quasilocal (Fig. [TT]) functional. 

Another fascinating collective phenomenon that can presently be described for the non-local 
[291 130] and quasilocal functional [31] is the fission of very heavy nuclei. In Fig. [T2J, an example 
of fission-path calculations performed in 258 Fm is shown in function of the elongation and 
shape- asymmetry parameters. One obtains correct description of the region of nuclei where the 
phenomenon of bimodal fission occurs and predicts regions of the trimodal fission, see Fig. [T3j 

In recent years, significant progress was achieved in determining the multipole giant resonances 
in deformed nuclei by using the RPA and QRPA methods. In light nuclei, the multipole modes 
can be determined for the nonlocal, relativistic, and quasilocal functionals, see Figs. [HI [T5l and 
[T6l respectively. In heavy nuclei, such calculations are very difficult, because the number of 
two-quasiparticle configurations that must be taken into account grows very fast with the size of 
the single-particle phase space. Nevertheless, the first calculation of this kind has already been 
reported for 172 Yb, see Fig. [T71 The future developments here will certainly rely on the newly 
developed iterative methods of solving the RPA and QRPA equations [3H [35]. 

The EDF methods were also recently applied within the full 3D dynamics based on the 
time-dependent mean- field approach. In Ref. [10], the spin-independent transition density was 
calculated in the 3D coordinate space for the time-dependent dipole oscillations. It turned out 
that one of the Steinwedel- Jensen's assumptions [H], 5p n (r;t) = —8p p {Y]t)^ was approximately 
satisfied for 8 Be. In contrast, in 14 Be, large deviation from this property was noticed. Figure [18] 
shows how transition densities 5p n (r;t) (lower panels) and 5p p (r;t) (upper panels) evolve in 
time in the x-z plane. The time difference from one panel to the next (from left to right) 




Figure 10. Energies and the average axial 
deformations for the two lowest quadrupole states 
with angular momenta + , 2 + , 4 + , and 6 + in 
24 Mg, together with the mean-field (dotted) and the 
corresponding angular-momentum-projected energy 
curves and squares of collective wave functions. 
Reprinted figure with permission from Ref. |27J. 
Copyright 2010 by the American Physical Society. 



Figure 11. Excitation spectra and 
B(E2) values in e 2 fm 4 compared to the 
available experimental data in 24 Mg. 
The spectrum is subdivided into a 
ground-state band, a 7 band, and 
additional low-lying states that do not 
necessarily form a band. From Ref. |28j; 
picture courtesy of M. Bender. 




Figure 12. The energy surface of Fm 
calculated for the quasilocal EDF as a function 
of two collective variables: the total quadrupole 
moment Q20 representing the elongation of 
nuclear shape, and the total octupole moment 
Qso representing the left-right shape asymmetry. 
Indicated are the two static fission valleys: 
asymmetric path aEF and symmetric-compact 
path sCF. From Ref. [32]; picture courtesy of A. 
Staszczak. 



Figure 13. Summary of fission pathway 
results obtained in Ref. [31]. Nuclei 
around 252 Cf are predicted to fission along 
the asymmetric path and those around 
262 No along the symmetric path. These 
two regions are separated by the bimodal 
symmetric fission around 258 Fm. In a 
number of the Rf, Sg, and Hs nuclei, all 
three fission modes are likely (trimodal 
fission). From Ref. [3T] 



roughly corresponds to the half oscillation period. White (black) regions indicate those of positive 
(negative) transition densities. One sees that significant portions of neutrons actually move in 




Figure 14. Dipole excitations 
24 Mg determined for the nonlo- 



m 



cal Gogny D1S functional. From 
Ref. |36j; picture courtesy of S. Peru. 



Figure 15. Dipole excitations in 20 Ne determined 
for the relativistic NL3 functional. Reprinted figure 
with permission from Ref. [37]. Copyright 2008 by the 
American Physical Society. 
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Figure 16. Quadrupole excitations in 24 Mg deter- 
mined for the quasilocal Skyrme SkM* functional. 
Reprinted figure with permission from Ref. [38] . 
Copyright 2008 by the American Physical Society. 
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Figure 17. Dipole excitations in 172 Yb 
determined for the quasilocal Skyrme 
SkM* functional. From Ref. [39]; 
picture courtesy of J. Terasaki. 



phase with protons. 

An interesting 3D EDF time-dependent calculation was recently performed for the a- 8 Be 
fusion reaction [42]. Although this calculation aimed at elucidating properties of the triple-a 
reaction, it was performed at the energy above the barrier, where the time-dependent mean-field 
approach can lead to fusion, whereas the real triple-o: reaction involves tunneling through the 
Coulomb barrier. Nevertheless, the studied tip-on initial configuration in the entrance channel, 
shown in the upper panel of Fig. [T9l is probably the preferred one as it must correspond to 
the lowest barrier. The calculations lead to the formation of a metastable linear chain state 
of three a-like clusters which subsequently made a transition to a lower-energy triangular a- 
like configuration before acquiring a more compact final shape, as shown in the lower panels of 
Fig. EH 




Figure 18. Nuclear TDHF transi- 
tion densities in 14 Be, see the text. 
Reprinted from Ref. [2D], Copyright 
2007, with permission from Elsevier. 




Figure 19. Nuclear TDHF densities in 
14 Be, see the text. Reprinted figure with 
permission from Ref. |42j. Copyright 2010 
by the American Physical Society. 



4. Conclusions 

The EDF methods are presently very intensely studied and developed in nuclear physics. Apart 
from the subjects covered in the present short review, there is a number of topics that could not 
be discussed here, such as: formal aspects of the DFT for self-bound systems; symmetry breaking 
effects; ab initio derivation of the EDF from the chiral perturbation and Brueckner-Hartree-Fock 
theories; studies of weakly-bound systems and the continuum phase-space effects; functionals 
non-local in time and the adiabatic connection; new-generation functionals with higher-order 
derivatives and/or richer density dependencies; self-interactions and self-pairing in the EDF; and 
ambiguities and inconsistencies when extending the EDF methods to multi-reference applications. 

Moreover, there are several aspects related to the methodology that were not covered here, 
such as: the EDF methods based on natural occupation numbers and/or orbitals; second, 
extended, and/or self-consistent (Q)RPA methods; proton-neutron interactions and isovector 
terms in the EDF; extrapolations to exotic nuclei and astrophysical applications; equation of 
state, symmetry energy, and neutron stars; relations between functionals describing infinite 
systems and finite nuclei; and adjustments of parameters, confidence intervals, and correlations. 

There is also a number of other interesting applications of the nuclear EDF methods, such 
as: description of tensor effects and the spin-orbit splitting; functionals describing pairing 
correlations; particle-vibration coupling, single-particle spectra, and widths of giant resonances; 
fusion barriers and cross sections within static and time-dependent calculations; fermion systems 
in the unitary regime; neutron skins and pygmy resonances; time-odd terms versus spin and 
orbital Ml resonances, spin-isospin resonances, and particle- vibration coupling and polarization; 
incompressibility, effective mass, and monopole resonances; cluster structures and models; 
chirality in rotational bands; di-neutron correlations and deformation in nuclear halos; and 
Coulomb frustration effects in the superheavy nuclei and crust of neutron stars. 

In general, the EDF methods provide us with universal understanding of global low-energy 
nuclear properties and feature an impressive array of applications. These methods can be 
rooted in the effective-theory approach whereupon the low-energy phenomena can be successfully 
modeled without resolving high-energy properties. Further progress strongly relies upon the use 
of high-power computing and faces the challenge of working out a consistent scheme of consecutive 
corrections that would allow for the increased precision and predictive power. 
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